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Abstract 

While static equilibria of flexible strings subject to various load types 
is well understood textbook material, the combinations of the very same 
loads can give rise to complex spatial behaviour. We show that flexible 
strings subject to combined smooth loads may not have smooth solutions 
in certain ranges of the load ratios. This nonsmooth phenomenon corre- 
sponds to the collapse geometry of inflated tents. After proving the nonex- 
istence of smooth solutions for a broad family of loadings we identify two 
alternative, critical geometries immediatly preceeding the collapse. We 
verify these analytical results by dynamical simulation of flexible chains 
as well as with simple table-top experiments with an inflated membrane. 
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1 Introduction 



The geometry and equilibrium configuration of flexible strings has been studied 
ever since the classical papers on the catenary by Bernoulli [2 , Leibniz [4j 
and Huygens ^ in 1691. In addition, Jacob Bernoulli discovered that the same 
curve is the solution the problem of the shape of a sail under Newtonian parallel 
wind [3 . Classical textbooks [7J extend the description of stable equilibria of 
flexible strings to various load types (e.g. gravity, wind or hydrostatic pressure) 
as analytical solutions of boundary value problems associated with ordinary 
differential equations. Such stable configurations can only fail due to intrinsic 
material nonlinearities [8 . On the other hand, elastic structures have additional 
failure modes due to global geometrical instabilities, as described first by Euler 
in 1744 [10 . More recently, the nonlinear behaviour of strings has been discussed 
by Antman [1 who describes three types of loading: gravity {g) , suspension 
bridge load {p) and hydrostatic pressure {h). In addition we also consider the 
parallel newtonian wind w discussed by Jacob Bernoulli and henceforth we refer 
to these as classical loads. Antman identifies in case of classical loads multiple, 
smooth solutions for nonlinearly elastic strings. One of these solutions is stable 
under tension and there exist multiple, smooth compressive solutions which 
are unstable. Here we look at the very same four classical loading cases and 
show that their linear combination, resulting in non- classical loads, can lead to 
unexpected and complex phenomena. While the governing ODEs can be readily 
derived, their solutions are nontrivial. In fact, for these non-classical loads we 
show that for certain finite ranges of the load ratios p/g^h/g^w/g no smooth 
solution exists. 



2 Mathematical background 
2.1 Governing equations 

We assume a flexible (ie. zero bending stiffness), inextensible string in a vertical 
plane, parametrized by the arclength s with slope a{s) measured with respect 
to the horizontal axis x, and tension T{s). We start with the formula for the 
hoop stress which, according to Gordon [13 originates in Mariotte's formula for 
cylindrical vessels: 

/« = I (1) 

describing the relationship beteen the radial force (normal pressure) /n, the 
tension T and the radius of curvature R. Equation ([T]) may be rewritten as 

n = a=^ (2) 

where denotes the curvature and () = d/ds. By supplementing the radial 
(normal) equilibrium equation Q with the (trivial) tangential equilibrium 

T = ft (3) 
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and noting that the fiexural rigidity did not enter into these equations, we see 
that ([2][3| provide the full description of the equilibrium of flexible strings. This 
seems to be of some interest as Mariotte's formula appears to be from 1680, 
whereas the equations for the string were first published (simulatneously) by 
Johann Bernoulli [2 , Leibniz [4| and Huygens [5 over a decade later. 

2.2 Cartesian frame 

If we are also interested in the spatial configuration x{s)^ y{s) in the global [xy] 
orthogonal frame then we have to supplement (2][3) by 



y = sin(a) (4) 
X = cos(a), (5) 

(expressing inextensibility) with corresponding boundary conditions 

x{0) = (6) 

^(0) = (7) 

x{L) = xo (8) 

y{L) = yo (9) 

which, together with the governing equations define a well posed boundary value 
problem (BVP) ([2][9|. Without loss of generality we assume xq > and for later 
reference and we introduce the 'boundary slope' as 

- ^ < ao = arctan ( — ) < ^. (10) 
2 \xo J 2 

Before continuing, we show that the connection to the Cartesian represen- 
tation is straightforward. Indeed, since 

cos(a) = — ==, = 77V"^2A372' Ql' = d/dx, (11) 

and the horizontal component H of the tangential force T is obtained as 

H = Tcos{a), (12) 
equations ([2|-(|3| may be rewritten as 

// fn /-I Q\ 

ii cos^(a) 

H' = ft-fnV'. (14) 

If we now consider a vertical load ^, so the normal pressure and tangential load 
fn , ft can be expressed as 

fn = gcos{a) ft = gsm{a). (15) 
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By substituting (15) into (13)-(14) we have 



y 



H' = 

from which the classical ODE for the catenary 



y 



H 



(16) 
(17) 

(18) 



immediately follows. If instead of the dead weight g we substitute the 'bridge 
load' p = gcos{a) then we get 



fn = pcos^{a) ft = pcos{a) sin(a) 



and by substituting into (13)-rtl4h we arrive at the othe well-known ODE 



y 

H' 



p_ 
H 




(19) 

(20) 
(21) 



yielding the parabola. The "Newtonian" model of vertical wind load with inten- 
sity w/2 assumes elastic collision between the incoming particles and the string, 
resulting in loading normal to the string's tangent: 



/n = ^cos^(a), A=0, 



so substitution into ( 12 )-( 13 ) yields 

y" - 

r = 



Tcos(a) 




(22) 

(23) 
(24) 



which shows complete analogy to ( 16 ) and also yield the classical ODE for the 
catenary 

(25) 



y = ^vi 
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Thus the physical shape in the [xy] plane of the string under gravity and vertical, 
Newtonian wind is identical, however, in the first case the horizontal component 
of the tension is constant while in the sacond case the magnitude of the tension 
itself, so in the [a,T] phase space of (2]|3) the solutions appear as different 
trajectories. In case of hydrostatic pressure h we have 



fn = h, ft = 0, 



so substitution into (13)-(14) yields 

y" = 

H' = 



H 



ii + y") 



-hy'. 



(26) 

(27) 
(28) 
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2.3 Critical solutions of the IVP and their interpretation 



In general, the load components fn^ ft may depend on the spatial location, 
however, here we only treat loads which only depend on the slope a, i.e. we 
have 

/„ = /„(«) ft = ft{a). (29) 



We remark that all four special loads (gravity ( |T5| ), bridge load (19), Newtonian 
wind (22) and hydrostatic pressure (26)) satisfy (29). This implies the system 
([2][3| becomes autonomous and can be integrated (at least numerically) as an 
initial value problem (IVP) to find a{s)^T{s). Also, we observe that the IVP is 
solved by the simple Ansatz 



a{s) 



where at are the real roots 



constant. 



1,2,. ..n 



fn{a) = 



(30) 



(31) 



and we call (30) critical solutions of ([2j[3|). These solutions appear as straight, 
horizontal lines in the [a,T] phase space of (|2]|3| (see Figure [T]). Since the phase 
space is a unique representation of solutions of the IVP (i.e. trajectories cannot 
intersect), this implies that if n > 1 then the phase space is partitioned into 
invariant domains (vertical stripes, defined by consecutive roots at^at_^i.) All 
non-critical, smooth solutions will be trapped in one of these stripes, i.e. their 
slopes will be bounded both from below and from above. Of special interest is 
the case when the boundary slope ao (defined in ([lO|) is trapped in one of these 
stripes, i.e. 

3k : at ^ ao < at_^i (32) 
If ( 32 ) is fulfilled and the width of the stripe satisfies 



(33) 



then we call the relevant roots of (31) and in the next subsection 

we show that the existence of relevant roots excludes smooth solutions for a 
sufficiently long string. 



2.4 Nonexistence of smooth BVP solutions 

In this subsection we show that in some cases the BVP may have only non- 
smooth or self-intersecting solutions. We start with a simple observation about 
the latter. Let r{s) = [x{s),y{s)]^ be a smooth curve, where s G [0,L] (L > 1) 
denotes the arclength and we have the boundaries r(0) = [0, 0]^, r(L) = [1,0]^. 
Then we have 

Lemma 2.1. If Js"^ G [0, L] such that r(s'^) = [1, 0]-^ then t{s) is self-intersecting. 

Proof. We start by noting that at any local extremum of y{s) we have t{s) = 
[±1,0]-^, and then sketch a proof of the Lemma by contradiction. The denial of 
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Figure 1: Solutions in [a,T] phase space. Straight, vertical lines represent crit- 
ical solutions cooresponding to the roots of (31). Regular solutions are trapped 
in the stripes between critical solutions. Relavant roots exist if two critical so- 
lutions a^, Q^^^i are bracketing the boundary slope ao and Sk = —(^l<7r. 



Lemma 2.1 is equivalent to the statement that all extrema of y{s) are associated 
with X = — 1, however, r{s) has no self-intersection. Let the y{s) assume its 
global minimum and maximum in 5 G [0,L] at si,52 respectively and without 
loss of generality we assume Si < S2. Then the segment s G [^1,52] of r(s) will 
separate the stripe y G [y{si)^y{s2)] into two disjoint domains (cf. Figure [2|. 
Depending on which of the boundaries r(0),r(L) is contained in which of the 
disjoint domains, either r(0) and r{si — e) or r{s2 + e) and r(I/), or both pairs 
are separated by the segment s G [^1,52] of r(s), i.e. the curve can only be 
connected smoothly by self-intersection. □ 

We proceed to formulate our main claim, and as a preparation we rotate the 
global [xy] frame by 

(^ = 0.5K+a^+i) (34) 
and we refer to the coordinates in the rotated system as ^, a etc., so we have 

= -^fc+i- (35) 

Now we can formulate 



Theorem 2.2. If (31) has relevant roots and L > xq/ cos(a^) then the BVP 



has no smooth solution free of self-intersection. 
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Figure 2: Illustration for the proof of Lemma |2.1[ Segment s G [si, 52] of t{s) be- 
tween the global minimum and maximum separates the stripe y G [y{si)^y{s2)] 
into two disjoint domains. Observe that either r(0) and r{si — e) (case a) or 
r(s2 + e) and r(L) (case b), or both pairs (case c) are separated the segment 
s e [51, 52] of r{s). 



Proof. We again prove the statement by contradiction, i.e. we start by as- 
suming that a smooth, non self-intersecting solution a{s) exists. The lack of 
self-intersection implies, via Lemma jSTT] that 



3s'^ such that a{s'^) = ao, (36) 

i.e. any non self-intersecting smooth solution a{s) is confined to 

a{s) e {al.al^^). (37) 

from which (due to the symmetrization and due to the restriction Sk < tt) 
it also follows that 

cos(a(5)) > cos{al) = cos(q;^^;l) > (38) 
We can now express the arclength: 

^ dx Xo 



L= ds= ^3v<^^. (39) 
Js=o Jx=o cos(a) cos(a^) 

Since this contradicts the condition of the Theorem, our proof by contradiction 
is completed: we showed that if (31) has relevant roots then the BVP has no 
smooth solutions without self-intersections if the length L is above the critical 
value indicated in the Theorem. □ 
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2.5 Existence of nonsmooth and self-intersecting solutions 

Theorem 12 . 2 1 shows the non-existence of smooth solution free of self-intersections 



if (31) has n > 1 real roots and the distance between the relevant roots is 
less then tt. Nevertheless, non-smooth solutions can, in principle, exist in this 
case. These solutions consist of segments joined by kinks, i.e. a continuous 
curve where the derivative is discontinuous at some places. These solutions are 
physically admissible only if the equilibrium conditions are satisfied at the kinks. 
This can be achieved in three ways (see Figure |3|: 

(a) If the tension T does not vanish at a kink, the kink represents a jump of 
TT in the tangent and string segments in tension are joined to string segments in 
compression. While this certainly satisfies continuity and equilibrium, these so- 
lutions are hardly relevant as they contain compressed string segments, resulting 
in infinite degrees of instability. 

(b) If 7^ then we have non-constant tension and the tension force T may 
vanish at isolated points. If the kink with a 7r-jump in the tangent is located 
at this point then both segments (before and after the kink) may be in tension. 
This situation is physically relevant and indeed we can observe such kinks in 
the simulation for hydrostatic loading combined with gravity. Needless to say, 
in such a case the string segments before and after the kink may come into 
contact which is not included in our equations (however, it will be included in 
the numkerical model). 

(c) Another possibility is the degenerate case when the tension T identically 
vanishes along the whole string, in this case kinks with arbitrary angles are 
admissible. This scenario may be realized in the case of the gravity-bridge load 
combination at special value of the load ratio. The numerical simulation shows 
indeed string shapes containing several kinks. 



We can conclude that if (31 ) has relevant roots then the BVP has no smooth 
solution and non-smooth solutions may only exist under special conditions. Next 
we show that while for all four classical loads all roots are irrelevant, pairing 
either the bridge load or the hydrostatioc pressure with gravity produces relevant 
roots and nontrivial physical behaviour. 

3 Analytical examples 

3.1 Classical loads: no relevant roots 

We describe four well-known load types: 

1. Gravity, with uniform weight g per unit arclength 

2. "Bridge load" with uniform vertical load p per unit horizontal projection 

3. Newtonian parallel, vertical wind w per unit of horizontal projection 

4. Hydrostatic pressure with uniform load h per unit arclength, normal to 
the string's tangent 
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S=So 



T-(So> r(So)=-T-(So) 




S=So 



T(So)=r(So)=0 




S=Sn, T1So) 




T(s) 



T(So)=P(So)=0 3 



T(s)=0 



S=Sn S 



Figure 3: Nonsmooth solutions: existence of kinks, a) non-phyiscal kink where 
tensile and compressed string are joined, b) physical kink at transversal zero of 
T{s) c) physical kink when T{s) = 0. 



and in each case we set yo = 0. As we already mentioned in the Introduction, 
three load types (g^p^h) have been described in detail by Antman [1 . The 
equations for all four mentioned load-types are very well known and can be 



readily derived from ([2|3|. By substituting any of ([15|,([19|,(|22|,([2^ into ([31) 
we find that in the interval [— tt, tt] for gravity, bridge and wind loads we get two 
roots, at distance 5k = tt and in case of the hydrostatic pressure there are no 



roots (see Table [T] lines 1-4.) According to the condition Sk > tt in Theorem 2.2 



the two roots at distance tt in case of the vertical loads do not impose a constraint 
on the length L, so we can expect smooth solutions. However, these roots can 
be physically interpreted by saying that under these loads the string is either 
everywhere vertical or it has no vertical tangent. This statement agrees with 
Proposition 2.25 (p592) of Antman 

3.2 Coupled loads: existence of relevant roots 

We now investigate coupled loads , i.e. when either p^w or h is acting opposite 
to the gravity g. For the combined loads we introduce 

±0p = arccos(^/p), \0p\ < tt (40) 
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i^^y = arccos(^/iL'), \0w\ < (41) 

±0h = arccos(V^), \0h\ < tt (42) 

and as shown in lines 5-7 of Table [l] pT] ) may have multiple roots in [— 7r,7r], 
with distance 5^ < tt. 







f — 






i 




g cos(a) = 


i 


— tt/z 








o 
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1 ^ /o 


o 
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27 — \ 

pcos [a) = 
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—tt/z 
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+7r/2 
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h = 






5 


g-p 


cos{a){g — pcos(a)) = 


1 










2 


-7t/2 








3 


—Op 








4 


+0p 








5 


+7r/2 








6 


+(7r - 0p) 


6 


g-w 


cos(a)(^ — wcos{a)) = 


1 


-(7r-^») 








2 


-tt/2 








3 


— 0w 








4 


+0W 








5 


+7r/2 








6 




7 


g-h 


^cos(a) — h = 


1 


-Oh 








2 


+0h 



Table 1: Critical solutions for the loads g^p^w^h and the combined loads {g — 
p),{g — w)^{g — h). Observe that the ordering of the 6 roots in case of the 
gravity-bridge and the gravity-wind load combinations (lines 5-6) depends on 
the absolute value of 6p and respectively. The table shows the ordering for 
Op < 7r/2, i.e. g/p > O^g/w > 0. Here the distance = 26p < tt, so the 
conditions of Theorem 12.21 are fulfilled. 



3.3 Corollaries for coupled loads 

Using the results in Table [l] we can formulate three corollaries for these special 
load types. We observe that the relevant roots are symmetric with respect to 
zero, so in the symmetrization (34) we have (f = 0. Based on Theorem 2.2 and 
Table [T] we now have 

Corollary 1 If the load is a linear combination of gravity g and "bridge load" 
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(cf.table [T] line 5) then the BVP ([2]|9| has no smooth solutions if 



L p 

^ > - > 1 

xo g 



Corollary 2 If the load is a linear combination of gravity g and "wind load" w 
(cf.table [T] line 6) then the BVP ([2][9| has no smooth solutions if 

L w 

^ > — > 1 

xo g 

Corollary 3 If the load is a linear combination of gravity g and hydrostatic 
pressure h (cf.table [l] line 7)then the BVP ([2][9| has no smooth solutions if 

^<^<1 
L g 

These corollaries can be obtained from Theorem |2.2| by substituting the values 
from Table [l] for the specific load combinations. We can observe that in all three 
cases as the load ratios g/p,g/w,h/g pass through zero, the corresponding roots 
±6>p, ±6>^, ±0h cease to be relevant because 6k > tt. 

One can describe a bifurcation problem in all three cases, and for the bridge 
and wind loads these coincide, so henceforth we do not treat the wind load case 
separately. The natural bifurcation parameters for the bridge load and hydro- 
static cases are p/g and h/g, respectively. In Figure [4] we show the bifurcation 



diagrams, plotting the parameter values versus the roots of (31). We also in- 
clude the corresponding physical configurations in the [xy] plane as well as the 
corresponding phase portraits of ([2][3| in the [a, T] space. 

3.4 Limiting geometries 

The corollaries of the previous subsection indicate the intervals of the load- 
ratios where smooth solutions do not exist. The first natural question is to 
establish the limiting geometries corresponding to equilibrium configurations at 
the endpoints of these parameter intervals. 

3.4.1 Bridge and wind coupled with gravity 

In case of the gravity-bridge load combination we can observe that at the lower 
end of the critical interval of the load ratio, disp ^ g from above, the two relavant 



roots of (31) approach simultanously each other and zero. Also, we note 
that solutions at ±7r/2 exist, corresponding to vertical string segments. The kink 
between the horizontal and vertical segments can be only of type (b) (cf. Figure 
|3]/b) since the vertical segment has tension due to gravity. Consequently, we 
expect the limiting geometry of the smooth string to be fiat with almost vertical 
ends. On the other end of the critical interval the two relevant roots 0^3 4 = 
in row 5 of table [l] of (31) approach arccos(xo/i^), so we expect the limiting 



geometry to be a symmetrical, upward pointing wedge; here the kink is of type 
(c). 
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While both the bifurcation diagram and the hmiting geometries are identical 
for the bridge and the wind- type loads, the nonsmooth geometry inside the 
critical range is different. In case of the bridge load coupled with gravity, the 
resultant of both loads is vertical so, if the string is a straght line defined by one 
of the relevant roots 0^3 4 = then identically vanishing tension T{s) = is 
possible, so any number of type (c) kinks is admissible. 

In case of wind load coupled with gravity, the situation is different: here the 
resultant of the wind is normal to the string while gravity is vertical, so tension 
will change in straight solutions monotonically. This implies that only one end 
of a straight segment can have zero tension, so at most one kink is admissible 
if the string consists of straight segments. 

3.4.2 Hydrostatic pressure combined with gravity 

If we combine hydrostatic pressure with gravity, the result at the lower end of 
the critical parameter interval {h/g = xq/L) is analogous to the upper end of the 
previous case: the two relevant roots 2 — (in i"ow 7 of table ITl) approach 
arccos(xo/i^) so here we can see a downward pointing wedge witn a type (c) 
kink. The upper end of the critical interval {h/g = 1) is less trivial: here ( [31] ) 
has only one root, corresponding to the horizontal line. Unlike the previous case, 
here the vertical lines are not solutions so the boundary conditions can not be 
satisfied by using straight lines corresponding to the relevant roots. From this 
we can conclude that among smooth solutions a self-intersection must occur 
for h/g > 1. Since self-intersection is phyiscally inadmissible, the string will 
develop self-contact. Vertical segments of self-contact are in equilibrium and 
can be joined by non-straight segments. Inside the critical range predicted by 
Corollary 3 we expect to see smooth segments (convex from above) separated 
by vertical regions in self-contact. 

The critical geometries and their order is illustrated in Figure |4j 

4 Numerical simulations 

In order to study the behaviour of the system in the regime that is in between the 
limiting solutions (gray zone in Figure |4|, we carried out dynamical simulations. 
We modelled the flexible string using 100 beads, each connected to its two 
neighbours by very stiff springs with unit natural length and a steep repulsive 
interaction between the beads at close range to prevent the string crossing itself. 
External forces (but no torques) acted on the particles, and an artificial velocity- 
dependent damping was introduced in order to drive the system to equilibrium 
if one existed. The first and last beads were kept fixed. The equations of motion 
were integrated with the velocity- Verlet scheme [12 using a time step that was 
infrequently adjusted to keep the maximum displacement at each step under 
control. During the simulations, the gravitational force was kept constant, and 
the angle-dependent forces (bridge load or hydrostatic pressure) was increased 
by constant increments from zero up to a large value, where the string resembled 
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IVP in Bifurcation diagram BVP in EVP in Bifurcation diagram IVP in 

[a, T] space [a*, h/g] space [x,y] space [x,y] space [a*, p/g] space [a, T] space 



p/g 

Dr\rn UJ 




Load: hydrostatic (h) + gravity (g) 



Load : bridge (p) + gravity (g) 



Figure 4: Limiting geometries of the Boundary Value Problems for the com- 
bined gravity-bridge and gravity-hydrostatic pressure loads are shown in the 
two central (3rd and 4th) columns. The corresponding trajectories for the Ini- 
tial Valine Problems are shown in the first and last column. The bifurcation 
diagrams of the critical solutions are in the 2nd and 5th columns, with the stable 
and unstable solutions represented by solid and dashed lines, respectively. 



its limiting shape, and slowly decreased back down to zero again. In order to 
establish the existence of equilibrium shapes, the angle-dependent forces were 
kept constant at each value until such time that the largest velocity component 
was below a threshold, or (allowing for the possibility that no equilibrium shape 
is stable) a certain time has elapsed since the last change in the force magnitude. 

Our results are illustrated in Figure |5l the arrangement is analogous to that 
of Figure |4] For the hydrostatic pressure/gravity combination, self-contact was 
observed in the critical regime of the load ratio, while in the bridge load/gravity 
combination, no equilibrium solutions were observed here: the string displayed 
complex time-dependent dynamics. 
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h/g = 3.5 

h/g= 1.85 
h/g= 1.34 
h/g= 1.0 



h/g = 0.86 

h/g = 0.81 
h/g = 0.68 

h/g = 0.00 




p/g = 2.84 

p/g= 1.69 
p/g= 1.40 

p/g=1.39 
p/g=1.37 
p/g=1.07 



J p/g=1.06 



p/g = 0.00 



Figure 5: Results of numerical simulations of the flexible string under combined 
load. The left two columns correspond to the hydrostatic pressure/gravity com- 
bination, with h/g increasing (decreasing) in the left (middle) column. The 
right hand column corresponds to the bridge load/gravity combination. In the 
critical regime of load ratios, nonsmooth and self-intersection equilibrium solu- 
tions are observed in the hydrostatic pressure/gravity case, while for the bridge 
load/gravity combination we did not find any equilibrium solutions. 



5 Experiments: gravity and hydrostatic pres- 
sure 

We also carried out a simple table-top experiment for the combined gravity- 
hydrostatic pressure load case. Our goal was to identify the limiting geometries 



described in subsection 3.4.2| and to compare the experiment to the numerical 



simulations presented in section [4j Our equipment was an elongated rectangular 
container with dimensions 110x25x35 cm on top of which we applied a rectan- 
gular membrane of stress- free dimensions 110x40cm with an orthogonal mesh 
drawn for better shape identification. We assumed that in the middle segment 
of the elongated rectangle the cross section of the membrane will approximate 
the shape of the planar string. The internal pressure was regulated with a bicy- 
cle pump and a valve. Figure [6] shows a series of photographs taken during a full 
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load cycle. The wedge shape at the lower critical point as well as the relatively 
flat geometry prior to collapse are visible. The two segments in self-contact, 
joined by a kink, also resemble the numerical simulations. 




Figure 6: A table-top experiments on the combination of hydrostatic pressure 
and gravity. The pressure is increaseing from left to right from zero for the first 
six images, which correspond to the states shown in the left column (i.e increas- 
ing h/g ratio) of Figure [s] The last image shows the near-critical geometry as 
the load ratio is decreased from its maximum value. 



6 Conclusons and outlook 
6.1 Summary of results 

We presented a theory about the global geometry of flexible strings when sub- 
jected to two opposing smooth loads. Under the assumptions that load intensi- 
ties depend only on the sting's slope we showed that the key to the geometric 



collapse are the roots of the critical equation (31) and we proved that the ex- 
istence of relevant roots (defined in equations (32), (|33|) is the necessary and 
sufficient condition for the nonexistence of smooth equilibrium configurations. 
We investigated four classical loads (gravity, bridge load, wind, hydrostatic pres- 
sure) and showed that if we couple any of the latter three with the former (grav- 
ity) then relevant roots appear and the geometry becomes nonsmooth. While 
the behaviour under any one of these loads has been known for centuries, this 
appears to be an interesting addition to the classical string theory. 

Beyond proving the nonexistence of smooth solutions we also described the 



non-smooth regime. In all three investigated cases the roots of (31) underwent 
a saddle-node bifurcation as the load ratio was varied. One end of the critical 
regime with no smooth solution is the bifurcation point, the other end is deter- 
mined by the boundary conditions (length of string and horizontal span). We 
determined the smooth limiting geometries at both (lower and upper) ends of 
the critical regime. In case of bridge and wind load we found that the lower 
limiting geometry is a flat, rectangular shape while the upper limiting geometry 
is an upward pointed, symmetrical wedge. In case of hydrostatic pressure the 
lower limiting geometry is a downward pointing wedge, while the upper limiting 
geometry can not be determinded as self-intersection occurs before the end of 
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the critical regime can be reached. Nevertheless, the midlle portion of the string 
is approaching the horizontal tangent. 

By using numerical simulations we verified not only the above findings but 
also determined the nonsmooth geometries inside the critical range of the load 
ratios. In case of hydrostatic pressure we ran a simple table-top experiment 
which showed fair qualitative agreement both with theory and the numerical 
simulations. 



6.2 Open questions and potential generalizations to mem- 
branes 

While our results apply to the geometry of membranes only under very special 
conditions, the analogy between the models is apparent. The equilibrium equa- 
tions of thin membranes show a close analogy to ([2|-([3| if no in-plane shear is 
admitted: 



fn 

Y 

ki i = l,2 (44) 



(43) 



where R = ^2 is the sum of the principal curvatures, refers to the partial 
derivatives of the tension T and ft^i denote tangential loads [15], [M]. The main 
difference between the 2D problem (strings) and the 3D problem (membranes) is 
that while strings are intristically fiat (i.e. they do not have intristic curvature), 
the intristic geometry of membarnes is nontrivial and can be described by the 
Gaussian curvature. 

In case of intristically fiat membranes and simple boundary conditions (e.g. 
a rectangular membrane supported along two parallel edges) our results directly 
apply and this was also confirmed by our table- top experiments. 



In case of intristically curved membranes, equations (43)- (44) have to be 
supplemented by the Gauss-Mainardi-Codazzi equations [15 . While the inte- 
grabilty of this coupled system has been investigated in [14 and is beyond the 
scope of the current paper, we believe that some qualitative features of our 
results still apply. In particular, we expect that in case of intristically curved 
membranes the geometry of collapse will also include non-smooth shapes. 



Our assumption on the slope-dependence of loads in ( 29 ) is not necessarily 
true, in particular, strings with non-uniform mass distribution will have gravity 
loading which depends explicitly on the arclength. We expect our results to 
be structurally stable and this is supported both by the numerical simulation 



and by the table-top experiments. If assumption (29) is not valid then one can 
not solve ([2|-([3| by the simple Ansatz (30). However, we expect that for loads 



slightly varying with arclength and, possibly with x and instead of critical 
solutions with constant slope one will find critical solution with bounded slopes. 
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6.3 Potential applications 

Beyond giving insight into the global geometry of strings, our result can be also 
used to predict the approaching collapse of inflated tents in which hydrostatic 
pressure is used to act against gravity. In this case our model indicates that 
before collapse the middle of the roof will become horizontal and the curvature 
will decrease. The recent collapse of the inflated Metrodome in Minneapolis 
([11]) shows that this appears to be the case indeed. Here the collapse was 
caused by snow accumulating on the tent. 

It might be of interest to note that if tents were pressurized by a mechanism 
producing the Newtonian wind-type load then, as the weight increased (or the 
pressure decreased) then the slope towards the middle of the tent would increase 
and this would prevent an incremental collapse due to the accumulation of snow. 
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